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[57] ABSTRACT 

The combination of trellis coding and MPSK signaling 
with asymmetry (nonuniform spacing) to the signal set 
is disclosed with regard to its suitability for a fading 
mobile satellite communication channel. For MPSK 
signaling, introducing nonuniformity in the phase spac- 
ing between signal points provides an improvement in 
performance over that achievable with trellis codes 
symmetric MPSK signaling, all this without increasing 
the average or peak power, or changing the bandwidth 
constraints imposed on the system. Block interleaving 
may be used to reduce error and pilot tone(s) may be 
used for improving the error correction performance of 
the trellis decoder in the presence of channel fading. 
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TRELLIS CODED MODULATION FOR 
TRANSMISSION OVER FADING MOBILE 
SATELLITE CHANNEL 

5 

ORIGIN OF INVENTION 

The invention described herein was made in the per- 
formance of work under a NASA contract, and is sub- 
ject to the provisions of Public Law 96-517 (35 USC 
202) in which the Contractor has elected not to retain 10 
title. 

This application is a continuation of application Ser. 
No. 929,876, filed 11/13/86, now abandoned. 

BACKGROUND OF THE INVENTION 15 

This invention relates to transmission of trellis coded 
modulation (TCM) of digital voice and data, and more 
particularly to TCM in combination with multiple 
phase-shift-keyed (MPSK) signaling, and with asymme- 
try (nonuniformity) in the phase shift of the MPSK 20 
signal set for reliable high quality communications in 
the presence of fading conditions, particularly when 
limitations of power and bandwidth are imposed. 

There is a growing need for reliable transmission of 
high quality voice and digital data in satellite-based land 25 
mobile communication systems. These systems, which 
will be part of an emerging all-digital network, are both 
power and bandwidth limited. To satisfy the bandwidth 
limitation, one can employ bandwidth efficient modula- 
tion techniques such as those that have been developed 30 
over the recent past for terrestrial microwave commu- 
nications systems. Examples of these are multiple phase- 
shift keyed (MPSK) signaling, quadrature amplitude 
modulation (QAM) and the various forms of continuous 
phase frequency modulation (CPM). When power is 35 
limited, forward error correction (FEC) coding is ordi- 
narily used. 

When limitation of power and bandwidth are both 
imposed simultaneously, as in the mobile satellite appli- 
cation, it is most often not possible to achieve the de- 40 
sired data rate of 4.8 or 9.6 kilobits per second (kbps) 
with either technique acting alone. Instead, what is 
required is the integration of bandwidth efficient modu- 
lation scheme with some form of FEC coding to exploit 
the best possible attributes of both. 45 

In the past, coding and modulation were treated as 
separate operations with regard to overall system de- 
sign. In particular, most earlier coded digital communi- 
cation systems independently optimized: (a) conven- 
tional (block or convolution^) coding with maximized 50 
minimum Hamming distance, and (b) conventional 
modulation with maximally separated signals. 

About a decade ago, using random coding bound 
arguments, it was shown that considerable performance 
improvement could be obtained by treating coding and 55 
modulation as a single entity. (J. L. Massey, “Coding 
and Modulation in Digital Communications,” Proc. 
1974 Int. Zurich Seminar on Digital Commun., Zurich, 
Switzerland, March 1974, pp. E2(l)-(4).) Many years 
later, this concept was formalized into a rigorous theory 60 
which showed that optimally designed rate n/(n-f-l) 
trellis codes suitably mapped (to maximize Euclidean 
distance) into conventional 2" + 1 -point signal sets can 
provide significant coding gain without bandwidth ex- 
pansion when compared with uncoded conventional 65 
2"-point signal sets. (G. Ungerboeck, “Channel Coding 
with Multilevel/Phase Signals,” IEEE Trans, on In- 
form. Theory, Vol. IT-28, No. 1, January 1982, pp. 


2 

55-67.) It is this work that has laid the foundation for 
the design and development of all power and band- 
width efficient digital modems found in practice today 
and those that are to come in the future. 

The most common application of such trellis coded 
modulation (TCM) techniques is in the new generation 
of modems being developed for the telephone channel. 
Indeed, the present state-of-the-art is a rate 6/7, 8-state 
trellis coded 128-point QAM which is capable of trans- 
mitting 14.4 kbps over good quality (Dl-conditioned or 
better) leased telephone lines. (J. Payton and S. Qureshi, 
“Trelhs Encoding: What it is and How it Affects Data 
Transmission,” Data Communications, May 1985, pp. 
143-152.) Thus if it is practical to send 14.4 kbps over 
the telephone channel, transmitting 4.8 or 9.6 kbps in- 
formation for a 5 kHz satellite channel (typical of pres- 
ent considerations) might appear to be simple. 

Several reasons make this supposition untrue most of 
which relate to the additional sources of degradation 
present on the mobile satellite channel but absent on the 
telephone channel. First, Doppler frequency shifts due 
to mobile vehicle motion can be a serious source of 
performance degradation if not compensated. Second, 
the fact that the 5 kHz mobile channel is actually a slot 
in a frequency-division multiple access (FDMA) sys- 
tem, brings on the problem of interference due to en- 
ergy spillover from adjacent channels. Third, the satel- 
lite channel is inherently a nonlinear one primarily due 
to the holding and positioning aid (HPA) in the trans- 
mitter. Thus, one must either employ constant envelope 
modulations and operate at full power or, if using non- 
constant envelope, but bandwidth efficient modulations 
such as QAM, then the HPA operating point must be 
backed off in power to produce an approximately linear 
channel. 

The most serious source of impairment that does not 
exist on the telephone channel is the combination of 
multipath fading and shadowing, i.e., for reliable perfor- 
mance, the system must combat short fades and recover 
quickly from long fades. Fading, which for mobile satel- 
lite channels can be assumed to be modelled by a Rician 
distribution, not only introduces an error floor into the 
system but also makes the problem of carrier recovery 
more difficult. Depending on the ratio of direct and 
specular (coherent component) to diffuse (noncoherent 
component) signal power, one might even be required 
to employ differentially coherent or noncoherent detec- 
tion techniques, thus sacrificing the power saving asso- 
ciated with coherent detection. Finally, even if the 
above sources of degradation were absent, the power 
limitation imposed by the mobile satellite channel 
would preclude transmission at the high data rates 
achievable on the telephone channel. Also, whatever 
technique is decided upon must be able to be imple- 
mented and installed in the vehicle with a minimum of 
cost and complexity, perhaps two orders of magnitude 
less than that associated with a telephone channel 
modem. 

SUMMARY OF THE INVENTION 

In accordance with the present invention, the combi- 
nation of trellis coding and MPSK signaling with the 
addition of asymmetry (nonuniform phase spacing) to 
the signal set provides a viable modulation and coding 
scheme for a fading mobile-satellite channel. For 
MPSK signaling, introducing nonuniformity into the 
phase spacing between signal points provides an im- 
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provement in performance over that achievable with 
trellis coding and symmetric MPSK signaling without 
increasing the average or peak power, or changing the 
bandwidth constraints imposed on the system. 

Whereas trellis coded modulation transmitted over an 
additive white Gaussian noise (AWGN) channel has 
been previously considered, the object of this invention 
is use of trellis coded MPSK in the fading environment 
of a mobile-satellite communication system. The object 
is achieved by using a combination of trellis coding and 
asymmetric MPSK signaling. The effect of the fading 
on the phase of the received signal is fully compensated, 
either by tracking it with some form of phase-locked 
loop, or preferably with pilot tone calibration tech- 
niques. Thus, the combination of trellis coding and 
asymmetric MPSK will correct for degradation due to 
the effect of fading on the amplitude of the received 
signal. The technique of interleaving and deinterleaving 
is employed to further combat the loss of data bits due 
to fading. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a functional block diagram of a trellis en- 
coded transmission system embodying the present in- 
vention, and FIG. la illustrates a basic analysis model 
for the system of FIG. 1. 

FIG. 2 is a diagram illustrating in graph A symmetric 
8-PSK and in graph B asymmetric (nonuniform) phase 
spacing of signal points in 8-PSK modulation. 

FIG. 3 is a diagram illustrating set partitioning of 
asymmetric 4-PSK modulation. 

FIG. 4 is a diagram illustrating a trellis for MPSK 
signal assignment of 4-PSK modulation. 

FIG. 5 illustrates a pair-state transition diagram for 
the trellis diagram of FIG. 4 

FIG. 6 is a graph of bit error probability performance 
vs. Eb/No for rate £ trellis coded QPSK in the presence 
of Rician fading; 2 states, K=10; ideal channel state 
information. 

FIG. 7 is a graph of bit error probability performance 
vs. E^/No for rate £ trellis coded QPSK in the presence 
of Rayleigh fading; 2 states; ideal channel state informa- 
tion. 

FIG. 8 is a graph of optimum asymmetry parameter 
vs. Eft/No for rate £ trellis coded QPSK in the presence 
of Rician fading. 

FIGS. 9 and 10 are graphs which illustrate the analo- 
gous results of FIGS. 7 and 8 for the case of no channel 
state, namely bit error probability performance vs. 
E^/No for rate £ trellis coded QPSK in the presence of 
Rician fading; 2 states, K=10, and in the presence of 
Rayleigh fading; 2 states. 

FIG. 11 illustrates the results of the simulation of §, 
16 state trellis coded 8-PSK modulation (TCM) over a 
Rician fading channel. 

FIG. 12 illustrates the effect of the doppler frequency 
on the results of the stimulation of §, 16 state trellis 
coded 8-PSK modulation over a Rician fading channel 
with channel state information (CSI) and the same inter- 
leaving parameters as in FIG. 11. 

FIG. 13 is a table which summarizes the results of 
FIGS. 11 and 12. 

FIG. 14 illustrates in a chart, performance of trellis 
coded modulation vs. QPSK modulation with dual tone 
calibration technique (TCT). 


DESCRIPTION OF PREFERRED 
EMBODIMENTS 

Referring now to FIG. 1 which shows a block dia- 
5 gram of the present invention, input bits representing 
data or digitally encoded voice are passed through a 
rate n/(n-f 1) trellis encoder 10 which is normally im- 
plemented with a combination of n shift registers (the 
memory of the encoding operation) and appropriate 
10 mod-2 adders (XOR gates). The encoder output sym- 
bols are then block interleaved in an interleaver 11 to 
break up burst errors caused by amplitude fades of dura- 
tion greater than one symbol time. While in practice the 
depth of interleaving is finite and chosen in relation to 
15 the maximum fade duration anticipated, for the purpose 
of theoretical analysis, the usual assumption is made 
herein of infinite interleaving depth. This assumption 
provides a memoryless channel for which well-known 
bit error probability bounding techniques can be ap- 
20 plied. However, the simulation results do reflect a finite 
interleaving depth. Thus, simulation results will be 
slightly pessimistic when compared with those derived 
from theory. The primary purpose of the theoretical 
analysis is to indicate through an example the trend of 
25 the performance behavior as a function of the various 
system parameters leaving the actual numerical perfor- 
mance to be predicted by the simulations. 

Groups of n+1 interleaved symbols are mapped with 
a read only memory (ROM) mapper 12 into the MPSK 
30 asymmetric signal set according to the set partitioning 
method. (G. Ungerboeck, “Channel Coding with Mul- 
tilevel/Phase Signals,” IEEE Trans, on Inform. The- 
ory, Vol. IT-28, No. 1, January 1982, pp. 55-67.) The 
in-phase and quadrature components of the mapped 
35 signal point are pulse shaped by a digital pulse shaper 13 
to limit intersymbol interference (ISI) and modulated 
onto quadrature carriers by a quadrature modulator 14 
for transmission over the channel. If pilot tone calibra- 
tion techniques are used to recover the faded carrier at 
40 the receiver, then the pilot tone (or tones) must be 
added to the data-modulated IF signal as indicated by 
an adder 15 before transmission on a RF carrier by a 
transmitter 16 over a fading channel 18. 

At the receiver 20, the RF carrier is received and the 
45 IF signal is recovered. The faded noise-corrupted in- 
phase and quadrature signal components of the IF sig- 
nal are: demodulated by a quadrature demodulator 21 
aided by the pilot tone(s) from an extractor 22; q-bit 
quantized by a quantizer 23 g for soft decision decoding; 
50 and then block deinterleaved by a deinterleaver 24. The 
metric chosen for the Viterbi algorithm 25 in a trellis 
decoder 26, which uses a metric computer 27, depends 
upon whether or not channel state information (CSI) is 
provided. (J. Hagenauer, “Viterbi Decoding of Convo- 
55 lutional Codes for Fading-and Burst-Channels,” 1980 
Int. Zurich Seminar on Digital Commun., Zurich, Swit- 
zerland, March 1980, pp. G2.1-G2.7.) 

As indicated in FIG. 1, a measure of CSI can be 
obtained from the power in the recovered pilot tone(s). 
60 Furthermore, the number of bits of quantization, p, in a 
P-bit quantizer 23^ for this operation can be much 
smaller than the number of bits of quantization q in the 
Q-bit quantizer 23 g since the accuracy of the CSI has 
only a secondary effect when compared with that of the 
65 soft decisions themselves. Finally, the tentative soft 
decisions from the Viterbi algorithm 25 are stored in a 
buffer 28 whose size is a design parameter. In particular, 
for the case of voice transmission, the total coding/de- 
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coding delay must be kept below about 60 ms so as not 
to be objectionable to the listener. Thus, for a given 
input bit rate, the decoder buffer and interleaving frame 
sizes must be limited so as to produce at most a 60 ms 
delay. For simplicity, an infinite buffer is assumed in the 
analysis, whereas the simulations reflect a finite buffer 
in accordance with the delay limitation. 

In order to better understand the use of the term 
“asymmetric MPSK” as applied to the present inven- 
tion, reference is made to FIG. 2 which shows in graph 
A a symmetric 8-PSK signal set uniformly spaced in 
phase. Graph B then shows an asymmetric 8-PSK signal 
set nonuniformly spaced in phase. It is this spacing that 
is referred to as asymmetric 8-PSK, even though there 
is symmetry but not uniform spacing in phase. 

Derivation Of The Pairwise Error Probability Bound 

The basic analysis model for the system of FIG. 1 is 
illustrated in FIG. la. The box 30 labelled encoder is 
actually the combination of the trellis encoder 10 (with 
binary input and output) and the function of the mapper 
12. This combination can be described as a trellis en- 
coder with binary input and M-ary output. (R. Calder- 
bank and J. E. Mazo, “A New Description of Trellis 
Codes,” IEEE Transactions on Information Theory, 
Vol. IT-30, No. 6, November 1984, pp. 784-791.) In our 
notation, M=2 n+1 . Also, if the input bit rate is denoted 
by R b, then the symbol rate input to the channel is R&/n. 
We denote a coded symbol sequence of length N by 

X=(xt,x 2 , . . . , x„) (l) 


6 

independent random variables (r.v.’s) and hence the 
channel probabilities satisfy 


Pn<£ \x,Z)= * p(yn\xn>Zn) 

N= 1 

and 

10 N W 

W(2) = It <fan) 
n— 1 


For any coded communication system, the decoder 
15 33 uses a metric of the form m(y,x;z) if side information 
is available and m(y,x) if it is not. Whatever metric is 
selected, it is desirable from the standpoint of simplify- 
ing the decoding process that it have an additive prop- 
erty, namely that the total metric for a sequence of 
20 symbols is the sum of the metrics for each channel input 
and output pair. In terms of equations 1, 4 and 6, this 
takes the form 


25 


m(# x; z) 


2 m(y n , x„; z„) 
n= 1 


The maximum-likelihood metric 


30 


m(x,x;z)=ln pMy|x,z) 


(9) 


(10a) 


where the kth element of x, namely xjt, represent the 
transmitted MPSK symbol at time k and is a nonlinear 
function of the state of the encoder s* and the n informa- 
tion bits, u k, at its input, viz., 

x*=f(s*,ujt) (2) 

The next state of the encoder s/t+i at the time k~b 1 is a 
nonlinear function of the present state s* and the n-bit 
input u* which is mathematically described by 

S*+l = g(SfoUA:) (3) 

Corresponding to x, the channel outputs the sequence 

y=(yi» y* • • ■ > y n ) (4) 

where the kth element yjt, representing the output at 
time k, is given by 

y*=pArX*-Hifc (5) 


when side information is available or 

m(y,x)=In pMy ]x) ( 10b > 

35 

when no side information is available satisfies the re- 
quirement in equation 9. This is easily seen by substitut- 
ing equation 7 in equation 10 and recalling that the 
logarithm of a product equals the sum of the logarithms. 

40 Thus, we shall use this metric in what follows. Also for 
simplicity, we shall use only the notation of equation 
10a, keeping in mind that the case where no side infor- 
mation is available is a special case where the metric 
m(y,x) does not depend on z. 

45 To find an upper bound on the average bit error 
probability performance of the Viterbi decoder, we 
must first find the pair-wise error probability P(x— >x) 
which represents the probability of choosing the coded 
sequence x— (xi, X 2 , . . . , x/v) instead of x=(xi, X 2 , . . . 

50 , Xiv). Since the decoder incorrectly decides x is the 
transmitted sequence when 


In equation 5, pk is a normalized (unit mean-squared 
value) random amplitude whose probability statistics 
depend on the fading channel model (to be discussed 
below) and n*is a sample of a zero-mean Gaussian noise 
process with variance cr 2 . The factor pk is the fading 
amplitude for use as side information, z n , in the trellis 
decoder under fading conditions. If this side informa- 
tion is available, then the corresponding side informa- 60 
tion sequence z will be denoted by 


m(y,x;z) = m(y,x;z) (11a) 

55 or, equivalently, 

N N (lib) 

2 m(y n , x„; z n ) ^ 2 m(y„, x„; z n ), 

n—l n = 1 

then, the pair-wise probability P (x->x) is given by 


z=(zi,z 2 • • • , z n ). (6) 

As stated hereinbefore, we shall assume adequate 65 
(theoretically infinite) interleaving and deinterleaving 
so that the coding channel is memoryless as indicated by 
box 32 in FIG. la. Under this assumption, the pjt’s are 




( N 

= Pr{ 2 

V*=i 


m(y, 


N 


, x n ; z„) ^ 2 m(y nt x n ; z n ) \ x 
n = 1 


*) 


( 12 ) 


Applying the Chemoff bound, equation 12 becomes 
equation 10 
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-continued 


( 21 ) 


(13) 

As -*-D = w E{exp(\[m(y n , £ n ; z n ) - m(y„, x„; z n )])|*} 
nerj 

where “E” is the statistical expectation operation, X is 
the Chemoff parameter to be optimized, and tj is the set 
of all n such that x^x*. To simplify equation 13 fur- 
ther, we must specify whether or not side information is 
available. 

a. Ideal Channel State Information (CSI) 

The assumption of ideal channel state information is 
tantamount to assuming that the side information r.v. z„ 
is equal to the fading amplitude p n . Thus, since n n is 
Gaussian distributed and in general, x„ and y n are com- 
plex quantities, then using equations 5 and 10a gives 

m(y rt ,x n ;z n ) =s — | y„ — p„x n | 2 (14) 

where for simplicity we shall ignore the icr 2 factor since 
it would later be absorbed in the Chemoff bound pa- 
rameter X anyway. Substituting equation 14 into equa- 
tion 13, and conditioning on £, a fading correction fac- 
tor obtained by monitoring the pilot signal amplitude, 
we get 


A 

d 2 (&x) = 2 p n 2 \x n ~ x„l 2 

ntTj 

represents the square of the weighted Euclidean dis- 
tance between the two symbol sequences x and x. 

Finally, the pairwise error probability upperbound is 
obtained by averaging equation 20 over the probability 
density function (p.d.f.) of £. 

b. No Channel State Information (CSI) 

When no channel state information is available, the 
metric of equation 10b becomes 

m(y„,x ff )= - [ y„-x„ j 2 (22) 

Substituting equation 5 into equation 22, then, analo- 
gous to equation 15, we get 

(23) 

P(x-*x\p) = 7T exp [— A j x n — x n \ 2 ~ 2X(p n — 
n €7} 

l)RfH>„ ( x n - *„)*}] X E{exp[-2\p„Re{n n (x n - x„)*}]} 

Again using equation 16, equation 23 simplifies to 


5 


10 


15 


20 


25 


A* — ilp) ^ exp (-Ap„ 2 |x* - 

“ rt€7J 


( 15 ) 


30 


Pix 


x\p) ^ 77- exp [— A|x„ - -X/ 7 1 2 (1 — 2\cr i y l ) — 
net) 


(24) 


x n \ 2 ]E{exp[-2\p n Re{n n (x n - x„)*}]> 


2A (p n — l)i?e{xfl(xn — x ra )*}] 


By representing the complex noise r.v. n* in terms of its 35 
real and imaginary parts where re{n*} and Im{n*} are 
uncorrelated, zero mean Gaussian r.v.’s each with vari- 
ance crjv 2 =cr 2 /2, it can be shown that 

E{exp [— * 2Ap rt Re{n rt (x rt — x n )*}]} — exp 

[2A 2 p n 2 o*iv 2 |x„~x„| 2 ] ( 1 6) 40 

Substituting equation 16 and equation 15 gives after 
some simplification 


f\x->x\p) = n exp(-Xp n 2 \x„ - x n \ 2 (l - 2Ao ’AT 2 )] 
net] 

Since <rN 2 is independent of n, optimizing equation 5Q 
17 over the Chemoff bound parameter X yields 


\ 


opt — 


1 

4 o-m 2 


(18) 


55 


which, when substituted in equation 17, produces the 
desired result, namely 


(19> 60 

Ai n DPn 2 \xn-xn\ 2 . D — exp (— 1 / 8 cta^) 
ne 7j 

Finally, equation 19 can be rewritten as 

/fe — i|g) S (20) 65 

where 


For constant envelope signal sets such as MPSK where 
| x | 2 = | x | 2 , equation 24 can be further simplified by 
noting that 

| x-x 1 2 =2Re{x(x-x)*} (25) 

Thus substituting equation 25 into equation 24, and 
renormalizing the Chemoff parameter (i.e., replace X by 
2 Xcr n 2 \ we get the desred result analogous to equation 
17, namely, 

P(x-^x|p,A)>D c2 (x,x|p,A) (26) 

where D is again given by equation 19 and 

(27) 

^feii^A) = 2 4\(p n - \)\x„ - x n \ 2 

^ ne Tj 

- 4A 2 Prt\x n — x n \ 2 — 4\ 2 d 2 (x,x) 
ner} 


Note that unlike equation 17, equation 26 cannot be 
optimized over X to yield a constant value for this pa- 
rameter. Thus, in this case, we must first average over 
the fading distribution. 

Derivation of the Bit Error Probability Bound 

To derive the upper bound on bit error probability 
from the pairwise error probability bound, we follow 
the transfer function approach taken in M. K. Simon 
and D. Divsalar, “Combined Trellis Coding with 
Asymmetric MPSK Modulation,” JPL Publication 
85-24 (MSAT-X Report No. 109), May 1, 1985. In 
particular, we first find the unconditioned pairwise 
error probability by averaging equation 20 or 26 over 
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the p.d.f. of p. (This averaging will be denoted by an 
overbar.) When this is done, the pairwise error proba- 
bility can be expressed in the form: 

P(X— (28) 
where for the ideal channel state information case 


A = 2 p n 2 \x n - x n \ 2 
nerj 

and for the no channel state information case 


A — 2 4X(p n — X)|x„ — x n \ 2 (29b) 

ne T) 

For uniformity of notation, we shall drop the depen- 15 
dence of the pairwise error probability on X for the no 
channel state information case. 

In dealing with upper bounds of this type, it is conve- 
nient to work with a pair-state transition diagram where 2Q 
the pair-state S* and pair-information symbol U& are 
defined as 


A 

Sk — (**» h) 


Uk =(«fo Ufc) 

where s& and u* are the estimates of the state of the 
decoder and the information symbol, respectively. We 
are in a correct pair-state when Sk—$k in an incorrect 
pair-state when sj^s*. Substituting equation 2 into 
equation 29, we can write the latter as 

A = 2 8 2 (S n , U„) (31) 

neri 

where 

8 2 (S„, Un) = Pn 2 L fan, u n ) - A^n> u n ) \ 2 (32a) 


where S 2 (S«, U n ) is given by equation 34. For the fading 
case, we simply replace 

D 62 (S„, U„) by D 52 (S n , Un) 

where again the overbar denotes averaging over the 
p.d.f. of p and S 2 (S„,U n ) is given by equation 32a or 32a, 
as appropriate, depending on the availability of channel 
state information (CSI). 

Finally for the case of no channel state information, 
we must minimize the upper bound of equation 35 over 
the Chemoff parameter to obtain the tightest upper 
bond. Recall that for the ideal channel state information 
case, we were able to perform this optimization at the 
pairwise error probability level. 

Characterization of the Fading Channel 

For mobile satellite communication, multipath fading 
produces a received signal with an amplitude which can 
be modelled by Rician statistics with parameter K rep- 
resenting the ratio of the power in the direct (line-of- 
sight) and specular component to that in the diffuse 
component. We assume that the phase distortion pro- 
duced by the fading is fully compensated for either by 
tracking it with some form of phase-locked loop or with 
pilot tone calibration techniques. (J. McGeehan and A. 
Bateman, “Phase-lock Transparent Tone-in -Band 
(TTIB): A New Spectrum Configuration Particularly 
Suited to the Transmission of Data over SSB Mobile 
Radio Networks,” IEEE Trans, on Commun., Vol. 
COM-32, pp. 81-87, January 1984). If shadowing is 
severe, then a Rayleigh statistical model becomes ap- 
propriate which can be looked upon as the limiting case 
of a Rician channel when K approaches zero. Of course, 
the case of no fading corresponds to a Rician channel 
with K approaching infinity. 

Mathematically speaking, the above statements cor- 
respond to a p.d.f. for the fading r.v. given by 


S 2 (S m U n ) = 4 \(p„ - \)\ A*n, u n ) - u„)\ 2 (32b) 

When no fading is present, equation 28 is given by 

p(x-4)=d a ( 33 ) 

where A is still of the form in equation 31 with, how- 


S 2 (S„,U„)= i f(s m u n )-f(s„,u„) 1 2 (34) 

In terms of the above definitions, it can be shown by 
analogy with the results in M. K. Simon and D. Divsa- 
lar, “Combined Trellis Coding with Asy metric MPSK 55 
Modulation,” cited above, that the average bit error 
probability P& is upper bounded by 

p^-jr <35) 6o 

where T(D,I) is the transfer function of the pair-state 
transition diagram whose branch label gains are modi- 
fied from those for the no fading case as follows. In the 
absence of fading each branch label gain has a factor 65 

D 52 (S m U„) 


2p(l + X)exp [-K - p 2 ( 1 + JO]/o(2p N K(l + K ) ) ; p 
0; otherwise 


where loOO is the zero-order modified Bessel function 
of the first kind. 

An Example 

Consider the case of rate \ trellis coded asymmetric 
QPSK using a 2-state trellis. The appropriate set parti- 
tioning is illustrated in FIG. 3, the trellis diagram in 
FIG. 4, and the pair-state transition diagram in FIG. 5. 
The performance of this system in the absence of fading 
was treated by M. K. Simon and D, Divsalar in “Com- 
bined Trellis Coding with Asymmetric MPSK Modula- 
tion,” cited above, with the following results: 


nm - 


• ; a » 4“ b = T’ d1 
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(37b) 


/£> 


Wfl = 

1 - ID 1+a 

where d as defined in equation 19 becomes 


and a is the ratio of powers between the I and Q chan- 
nels which is related to the angle <f> that defines the 15 
asymmetry (see FIG. 3) by 


. 1 T 

a — tan z — 


(39) 


Substituting equation 37 into equation 35 gives 


Pb ^ 


) l+a 


(40) 


<_Y 


Optimizing equation 40 over the assymetry produces 

Eb/No . (41) 


4 In D 

a = = ■ 


In 3 


, _ i 


which when substituted into equation 40 gives the de- 
sired upper bound 


(42) 


DP 2 P evaluates to - (43a) 

( Ui* Y, ij&ih' = _£L 

** ~ \ 1 + K + 07 ) U ’ 7 Wo 

which for the Rayleigh special case (K=0) becomes 


dp 2 & = • 


1 


I + 07 


60 


(43b) 


Evaluating T(D,I) of equation 37a using equation 43 
and performing the differentiation required in equation 
35 gives the upper bound on as 


Pb: 


, v (38) 

( E h \ E h ^ . 1° „ 

D as exp I — I; /V = system bit energy-to-noise fa = 

X °S o spectral density ratio 


Pi* 


: l + k + pry 




(i - fc r * 2 ) 2 

where 

1+A_ 


(44) 


(45) 


1 + K + Pry 


4 a 


/- i. 2. 3 ^1 - 4s 


To obtain the best performance in the presence of 
fading, one should optimize equation 44 over the asym- 
metry parameter a. Before doing this, however, we 
shall first examine the behavior of equation 44 for the 
20 symmetric case, i.e., a= 1, and the optimum asymmetry 
in the absence of fading as given by equation 41. 

Substituting a= 1 in equation 45, the parameters £,* £/ 
simplify to 


25 


= ■ 


4 K 


£2 = ?3 - 


1 + K + 4y 
IK 


;«i =■ 


! +* 


1 + * + 2y 


1 + /T + 4y 

l±K 


(46) 


fe“' 


1 + K + 2y 


A -J- exp(— 2Bj/iVo) 

a. Ideal Channel State Information 

Recalling equation 29a, we see that the transfer func- 
tion T(D,I) of the pair-state transition diagram for the 
case of ideal channel state information is obtained by 
merely replacing by 

DP 25 


in the branch label gains of equation 37a where /3=4, 5Q 
4/(1 +a), or 4 a/(l+a) as appropriate. For the Rician 
p.d.f. of equation 36, 


30 The curve labelled “symmetric” on FIG. 6 is a plot of 
the upper bound of equation 44 combined with equation 
46 as a function of the average bit energy-to-noise ratio 
E^/No with a Rician parameter K=10 (typical of the 
mobile satellite channel). 

35 When the value of a in equation 41 is substituted in 
equation 45, then after some simplification, equations 44 
and 45 can be written as 


40 


45 


fr'ft'expl-gf 4- ii)\ 
(1 - ?2'exp(-S2')) 2 


Pb ^ 

and 

^ ” 1 + K + P/ 
i = 1, 2, 3 

E b 
No 


B/K 




l±K 


(47) 


(48) 


I + K + Pi 


E b 

No 


55 


The behavior of equation 47 combined with equation 48 
is also illustrated in FIG. 6 by the curve labelled “opti- 
mum asymmetry (no fading).” 

For the Rayleigh case, the above results simplify 
even further. In particular, for the symmetric signal 
QPSK constellation, we get 


^ ( 1 + EbONo y \( _\ Y 1 ^ 

n * V Eb/2N 0 ) ^ 1 + Eb/No )\ 1 + tb/INo) ) 


(49) 


whereas for the optimum asymmetry in the absence of 
fading, 


Note that Eb represents the total average received 65 
power (direct plus specular plus diffuse) in the data 
bandwidth. In the general case of an n/(n+ 1) code rate, 

Efr in the definition of X would be replaced by nE^. 


P b ^ ^ V )( 1 + Eb/No )( 1 + Eb/No ~ In3 ) 


( 50 ) 
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These results are illustrated in FIG. 7. Note that here 
the curve labelled “optimum asymmetry (no fading)” 
gives worse performance than that of the symmetric 
case. Thus, at least here, we clearly see the need for 
performing the asymmetry optimization in the presence 5 
of the fading. 

To determine the optimum value of a for the Rician 
case, we need to differentiate equation 44 with respect 
to a and equate the result to zero. This leads to a tran- 
scendental equation which must be solved numerically. io 
Rather than doing that, it is more expedient to directly 
minimize equation 44 with respect to a using numerical 
techniques. When this is done, we obtain the optimum 
bit error probability bound labelled “optimum asymme- 
try (fading)” in FIG. 6. For K= 10, this curve lies quite 15 
close to the “optimum asymmetry (no fading)” curve. 
However, as we can already deduce form FIG. 7, this 
statement is not true for small values of K, in particular 
the Rayleigh channel with K=0. To exhibit the sensi- 
tivity of the optimum asymmetry condition to K, FIG. 20 
8 illustrates the optimum value of a as a function of 
E&/N0 with K as a parameter. 

For the Rayleigh case, we can indeed determine the 
optimum asymmetry condition in closed form. In par- 
ticular, differentiating equation 44 with K=0 in equa- 25 
tion 45 and equating the result to .zero has the solution 


-4 + Et/No (N 17 + 8 Eb/N 0 - 1) 
4(1 + Mb/m 


(51) 

30 


which when substituted back in equation 44 gives 


vi - < 52 > 35 

(3 + N 17 + SEb/No) [4 (Eb/No) z + 

2 

_ ^ (Eb/N 0 ) 2 (7 + Nl7 + 8 E b /No)] 

Ph = 11 

16(N 17 + ZEb/No - 1) (Eb/No) 2 ( 1 + Eb/No) 4 

40 

This result is illustrated by the curve labelled “optimum 
asymmetry (fading)” in FIG. 7 and is clearly superior to 
that corresponding to the symmetric signal constella- 
tion. 

Finally, for purpose of comparison, the correspond- 45 
ing upper bound on the performance of uncoded BPSK 
(same bandwidth as rate £ trellis coded QPSK) in the 
presence of Rician and Rayleigh fading is also illus- 
trated in FIGS. 6 and 7. The analytical results for these 
curves are well known and are given by 50 


l+JC ( KEt/No \ 

Pb = 1 + K + Eb/No exp ^ 1 + K + Eb/No ) 

for the Rician channel and 


(53) 


55 


result over this domain. In particular, for large E^/No, 
equations 49, 50, 52 and 54 respectively, become 


Pb* 


l 

(N 2 Eb /No) 2 


(symmetric) 


(49') 


PbZ 


(50') 


Y ln3 

> E b " 

A 1 + ln3 

) N ° . 


(optimum asymmetry — no fading) 


Pb ~ r- (optimum asymmetry — fading) 

(Eb/No) 2 


(52') 


and 


Pb ~ Eb/No ( uncoded) 54 j 

Thus, for example, comparing equation 49' with equa- 
tion 54', we see that the effect of coding is to change the 
rate of descent of the error probability vs. E&/No perfor- 
mance form an inverse linear to an inverse square law 
behavior. If the QPSK constellation is now designed 
according to the optimum asymmetry for no fading, 
then the performance is worse than that of the symmet- 
ric constellation by a factor of (1 -bln 3)/V2 In 3) or 1.3 
dB„ On the other hand, if the constellation is designed 
with the optimum asymmetry determined in the pres- 
ence of asymmetry, then, relative to the symmetric 
design, the performance is improved by a factor of 
square root of 2 or 1.5 dB. From FIG. 7, we see that 
these asymptotic results are almost achieved at an error 
rate of 10“ 5 . 

b. No Channel State Information 

Recalling Equation 29 b, we see that the transfer func- 
tion T(D,I) of the pair-state transition diagram for the 
case of no channel state information is obtained by re- 
placing DP this time by in the branch label gains of 
equation 37a where again ^3=4, 4/(1 + a), or 4a(/(l -ha) 
as appropriate. Unfortunately, for the Rician distribu- 
tion, the factor cannot be evaluated in closed form. It 
can, however, be expressed as a single integral with 
finite limits as follows: 


Wp = q~ k 


1 



T7(0)exp[7j 2 (0)]er/ci7(0)d0 


(55) 




Nl + K 


Pb ^ 


1 

1 + Eb/No 


( 54 ) This integral is easily evaluated using Gauss-Chebyshev 
techniques, i.e., 

60 


for the Rayleigh channel. 

_ Depending on the shape of the bit error probability vs 
E^/No curve, one can often deduce some important 
practical conclusions by examining the asymptotic be- 
havior of the curve. Since for the Rayleigh case (FIG. 65 
7), the error probability performance curves are essen- 
tially linear over a wide range of practical SNR/s, one 
can approximately apply the asymptotic (large E*/No) 


7 T (56) 

T](6)exp[7} 2 ($)]erfcri(0)dd sa 


7^ j y\{ek)^p[r\ 2 {Qk)]erfa]{e k ) 

where 0*4 (2k- l)ir/2N. 
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For the Rayleigh case, we can obtain a closed form 
result for this factor since, for K=0, n(0) becomes inde- 
pendent of 0. Thus 

D*Wp = 1 - (57) 

Evaluating T(D,I) of equation 37 a using equation 55 
and performing the differentiation required in equation 
35 gives the upper bound 


16X2(1 +2a)/(I+a) 

p/j :S min min 

X ^ 0 a <1 - £ 2 /)-16X2/(l + a))2 


a =r r- exp(4\ 2 Ei/N 0 ) 

MlkEtJ Nq) 2 


■ 2 

( ikEyNa \ 

A l +“ ) 

i 

( 2aX£t/N 0 \ 

A 1 + a ) 


( 4\ 2 Eb/Np ^ 
\ 1 + « ) 

( 4a\ 2 Eb/N 0 \ 

A J 


Substituting equation 63 into the transfer function of 
equation 37 a and performing the differentiation re- 
quired in equation 35 gives the approximate upper 
bound on P b (valid for large E&/No) 


it - e - "* £ 1 - ~^=r V&)exp(y t 2 (0))erfcTiA0)d0 J 


Pb ~ min - 
X 




_ mBb/2N 0 ) \nr , 

t?/< 0) = — vi -i.:aaa= cos0 

Nl + K 

and fin 7 j/= 1,2,3 are defined in equation 45. For the 3Q 
Rayleigh case one merely replaces \t and 17 ,<0) of equa- 
tion 59 by 

Si = 1 “ ^ Vi ex P(Vi 2 ) er f cr ti 35 

Vi = W0b/2N O ) 

and performs the same minimizations required in equa- 
tion 58. 40 

FIGS. 9 and 10 illustrate the analogous results to 
FIGS. 7 and 8 for the case where no channel state infor- 
mation is available. For simplicity of presentation, we 
have chosen not to illustrate the results for the values of 
optimum asymmetry determined from the no fading 45 
analysis since we have already made the point that 
asymmetry should be optimized in the fading environ- 
ment. Clearly, the lack of channel state information 
produces a noticeable degradation in system perfor- 
mance. To quantitatively assess this additional degrada- 50 
tion (at least for the Rayleigh channel), we now derive 
asymptotic results analogous to equations 49', 50', and 
52' for the no channel state information case. In particu- 
lar, we use the symptotic (large argument) expansion 55 
for erfc x, that is to say, 


. ^ 2?E b /N 0 y 


’ (t+T y 2E b/No ) 


Performing the minimization over X required in equa- 
tion 64 (actually the numerator is minimized since the 
denominator has little effect on this operation) gives 


{-Vrt-) 


which when substituted in equation 64 yields 


(-Hr 32 -) 


16{£y7V 0 ) 2 


Finally, the desired asymptotic results are 


KtEb/No? 


(symmetric) 


lEb/Np - M3 
Eb/Np — IriS 


■A 

/ l6(E b 


e 2 

WEb/Np) 1 


(optimum asymmetry — 


no fading) 


[1 -h x~ 2 ] 


P b < _ r- (optimum asymmetry — fading) 

MEb/Npy 


in which case equation 57 simplifies to 


2[X/3(£*/lV 0 )] 2 

Using the appropriate values of fi in equation 62, the 
branch gains of FIG. 5 become 


Simulation Results 

The results of a software simulation of the system 
block diagram of FIG. 1 which may be readily imple- 
65 mented by those skilled in the hardware art will now be 
described and presented. For the application at hand, 
the development of a simulation has manyfold purpose. 
First, it can be used to “verify” the theoretical results 
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set forth above, keeping in mind that the simulation is 
indicative of the exact system performance whereas the 
theoretical bit error rate expressions are upper bounds. 
Second, when the number of states in the trellis diagram 
becomes large (e.g., 16), determining the state transition 5 
diagram and its associated transfer function is a tedious 
task; in such cases, simulation is the more expedient 
approach. Finally, system degradation due to the finite 
size of interleaving and decoder buffer imposed by the 
practical constraint on the allowable total delay is ana- 
lytically intractable, particularly when coupled with 
that caused by the “noisy” carrier demodulation refer- 
ence produced by the pilot tone extractor. Hence, to 
predict true system performance corresponding to the 15 
real world environment, one must again turn to simula- 
tion. In the next paragraph, we expand upon the last of 
these issues. 

The block interleaver 11 of FIG. 1 can be regarded as 
a buffer with d rows which repressent the depth of 20 
interleaving and s columns which represent the span of 
interleaving. Thus, the size of the interleaver (in sym- 
bols) is d X s. Data is written into the buffer in successive 
rows and read out of the buffer (the order in which it is 
transmitted over the channel) in columns. At the re- 25 
ceiver, the block deinterleaver 24 performs the reverse 
operation, i.e., the received soft quantized symbols are 
written into the buffer in successive columns and read 
out in rows. In practice, the interleaving depth should ^ 
be chosen on the order of the maximum fade depth J 
anticipated which, for the fading mobile satellite chan- 
nel under investigation, depends on the doppler fre- 
quency or, equivalently, the vehicle speed. The smaller 
the doppler frequency, the longer the fade duration and 35 
vice versa. The interleaving span should be chosen on 
the order of the decoder buffer size. When this is done, 
the performance degradation (relative to that for the 
analytically tractable assumption of infinite interleaving 
depth and buffer size) will be inversely proportional to 40 
the product of interleaving size and doppler frequency. 

On the other hand, when pilot tone(s) are used for 
coherent demodulation as suggested in FIG. 1, then 
performance will degrade directly proportional to dop- 
pler frequency. The reason for this is that the bandpass 45 
filter(s) used in the pilot tone extractor 22 to isolate the 
piolet tone(s) from the modulation must have band- 
width sufficiently wide to include the doppler shift. 
Thus, the larger the doppler, the wider the bandwidth 
of the filter(s) and hence the “noisier” the extracted 50 
demodulation reference. Assuming infinite interleaving 
and decoder buffer size, one can use the same analytical 
approach as previously discussed to derive upper 
bounds on the bit error probability in the presence of 55 
the noisy carrier reference. In particular, we first find 
the probability density function (p.d.f.) of the phase 
error associated with the demodulation reference signal 
which itself depends on the envelope of the fading sam- 
ple. Next, the bit energy-to-noise ratio in the branch 50 
label gains (see FIG. 6 for example) should be degraded 
by the cosine of the phase error p.d.f., and then, as 
before, over the fading envelope p.d.f. Finally, using the 
averaged branch label gains in the transfer function 
bound T(D,I) gives the desired upper bound on average 65 
bit error rate. Even for the assumption of infinite inter- 
leaving and decoder buffer size, this computation is 
quite tedious. Thus, a simulation is preferable. 
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EXAMPLE 1 

This example is for the verification of the analytical 
results obtained in the previous section. In particular, 
the system block diagram of FIG. 1 was simulated for a 
rate J, 2-state trellis code, symmetric or optimum asym- 
metric QPSK modulation, and either a Rician or Ray- 
leigh fading channel. The interleaving size was chosen 
equal to 512 QPSK symbols (or equivalently, 512 input 
bits) which, for all practical purposes, approximates 
infinite interleaving. A doppler frequency of 100 Hz 
was chosen which makes the channel rapidly varying 
enough for the fading to be assumed independent form 
symbol to symbol. These two assumptions provide a 
memoryless channel as assumed in the analysis. Further- 
more, the buffer size was chosen equal to 32 bits which 
approximates the assumptions of an infinite bit buffer. 

Numerical simulation results for this example are 
superimposed as dashed lines on the analytical results of 
FIG’S 6 and 7 corresponding respectively to the cases 
of Rician and Rayleigh fading with ideal channel state 
information. The discrepancy between the solid and 
dashed curves reflects the looseness of the upper 
bounds but the relative behavior of the analytical curves 
compares well with that of the simulation results. 

EXAMPLE 2 

Here we consider the more practical case of a rate §, 
16-state trellis code combined with symmetric BPSK 
modulation. It was shown by M. K. Simon and D. Div- 
salar in “Combined Trellis Coding with Asymmetric 
MPSK Modulation,” cited above, that for this case the 
additional coding gain produced by the addition of 
asymmetry to the modulation is small and thus we have 
chosen to ignore it. Although, with much computa- 
tional effort and the assumption of infinite interleaving 
and buffer size, this system can be analyzed by the ap- 
proach taken in the previous section, our interest here 
lies in computing the performance with limited inter- 
leaving and decoder buffer size as follows. 

At the present time, this system is a candidate for 
National Aeronautics and Space Administration’s Mo- 
bile Satellite Experiment (MSAT-X) project whose 
objective it is to transmit 4800-9600 bps of digitally 
encoded speech over a 5 kHz RF channel with a bit 
error rate of 10“ 3 . To satisfy the previously mentioned 
constraint on total allowable delay, the interleaving 
size, interleaving depth, and decoder buffer size have 
been optimized at this bit error rate to achieve the mini- 
mum bit signal-to-noise ratio. For the specified delay 
constraint (60 ms), the size of the block interleaver and 
deinterleaver have been chosen equal to 128 BPSK 
symbols (or 256 input bits). With the above chosen 
interleaving size, the interleaving depth has been opti- 
mized by computer simulation and found to be equal to 
16 symbols. Thus, the interleaving span is 128/16=8 
symbols over the range of doppler frequencies from 20 
Hz to 100 Hz (vehicle speeds of 15 mph to 75 mph at 
UHF). Note, however, that for MSAT-X channels op- 
erating at low doppler frequencies such as 20 Hz, we 
can indeed have fade durations much longer than 16 
symbols. In this case, an interleaving size of 128 symbols 
is not sufficient and, as we shall see shortly, a significant 
performance penalty occurs. Finally, with the above 
delay constraint imposed, the buffer size was optimized 
through simulation and found to be 32 symbols (or 64 
bits). 
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FIGS. 11 and 12 illustrate the results of the simulation 
for perfect carrier and time synchronization and no 
intersymbol interferences (ISI). In particular, FIG. 11 
assumes a fixed doppler frequency of 100 Hz, a fixed 
Rician parameter K= 10, and various interleaving and 5 
channel state information options. Also shown as a 
reference point, is the performance of uncoded QPSK 
(same bandwidth as rate § trellis coded 8 PSK) at the 
chosen error rate of 10“ 3 . FIG. 12 shows the effect of 
doppler frequency on system performance for the case 10 
of ideal channel state information and the same inter- 
leaving parameters as in FIG. 11 . 

FIG. 13 summarizes in a table the above results by 
tabulating the required E^/No at a bit error rate of 10“ 3 
for each of the cases and also for Rician fading parame- 15 
ters of 5 and 6 dB. From these numerical results, one 
can assess the coding gain (reduction in required E&/N 0 
relative to uncoded QPSK) achieved in each case. For 
example, without interleaving and channel state infor- 
mation, transmitting § coded 8 PSK over the K=10 20 
Rician channel produces a coding gain of 1.6 dB. When 
128 symbol interleaving is added, the coding gain is 
increased to 3.1 dB. If, in addition to interleaving, we 
provide ideal channel state information, another 0.5 dB 
coding gain is achieved which brings the total coding 25 
gain to 3.6 dB in this particular environment. As is 
obvious from FIG. 11 , this coding gain would be 
greater at lower bit error rates. When the doppler is 
decreased from 100 Hz to 20 Hz, FIG. 12 and the table 
of FIG. 13 show a reduction in coding gain of 1 dB due 30 
to the limitation on the size of the interleaver. 

As previously mentioned, all of the results in FIGS. 

11 and 12 , and the table of FIG. 13 are for the case of 
perfect carrier synchronization. Using dual tone calibra- 
tion technique (DTCT), wherein two tones of equal 
power are inserted symmetrically at the edges of the 
data spectrum for the purpose of coherent demodula- 
tion, we have found by simulation that the noisy carrier 
reference produced by the appropriate pilot tone ex- 
tractor produces about a 2 to 3 dB degradation in per- 
formance depending on the value of the Rician parame- 
ter K. This is caused by the fact that the bandwidth of 
the pilots tone bandpass filters has to be chosen wide 
enough to accommodate the maximum doppler of 100 
Hz and the total power of the two pilots is 7 dB below 
that of the data (this ratio has been shown to be opti- 
mum). Simulation results with the DTCT technique are 
shown in FIG. 14 for Rayleigh and various Rician chan- 
nels under the assumption of ideal channel state infor- 
mation, 128 symbol interleaving, and a doppler fre- 
quency of 20 Hz (worst case for our application). Note 
that in this simulation, the bandwidth of the pilot band- 
pass filters was chosen to accommodate the 100 Hz 
doppler as would be the case in the actual hardware 
design. 

From the foregoing, it is evident that there is pro- 
vided a new and improved method of communication 
over a fading mobile channel which uses conventional 
trellis coding and conventional multiple phase shift 
keyed (MPSK) signaling with asymmetry (nonuniform) 
spacing in the phases of the MPSK signal set. This basic 
method may be readily implemented with conventional 
equipments using state of the art techniques. For exam- 
ple, the asymmetry in MPSK may be readily imple- 
mented using a read only memory (ROM) mapper. 
Trellis coded modems are known to tolerate more than 
twice as much noise as conventional modems at the 
same block error rate, and provide significantly higher 
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throughput by reducing the error rate by about three 
orders of magnitude, but that is through a nonfading 
telephone channel. The present invention provides reli- 
able, high quality digital voice and data transmission in 
a fading communication channel, such as a mobile satel- 
lite channel, where land mobile receivers are expected 
to experience severe fading conditions while in use as 
the mobile receiver moves through varying terrain. 

What is claimed is: 

1. A method of communication with a mobile re- 
ceiver over a fading mobile-satellite channel for digital 
coded voice and data comprising the steps of trellis 
coding said data at a transmitter, block interleaving said 
trellis coded data, and multiple phase-shift-keyed signal- 
ing the trellis coded and block interleaved data with 
constant amplitude and distinct phase spacing in a signal 
set of phases greater than four, and thereafter quadra- 
ture modulating the result on an IF carrier and modulat- 
ing the modulated IF carrier on an RF carrier for trans- 
mission through said mobile-satellite channel and in- 
cluding the step of inserting two pilot tones in the mod- 
ulated IF carrier for use in synchronous demodulation 
of the IF carrier at said mobile receiver, and further 
including the steps at said mobile receiver of extracting 
said pilot tones, using said pilot tones for synchronous 
demodulation of said IF carrier to coherently detect 
received trellis coded data, block deinterleaving de- 
tected trellis coded data, and decoding the deinter- 
leaved trellis coded data, said demodulating step includ- 
ing the steps of recovering information from extracted 
pilot tones for use as side information in said deinter- 
leaved trellis decoding step. 

2. A method as defined in claim 1 wherein said signal 
set is comprised of eight distinct phases. 

3. A method as defined in claim 1 wherein channel 
state information is derived from said pilot tones at said 
receiver in the form of a fading channel parameter, p n , 
for use as side information, z nt in the step of trellis de- 
coding, and the decoding process uses a metric of the 
form m(y„, x„; z n )= - |y«-p«x n | 2 for improving the 
error correction performance of the trellis decoding in 
the presence of channel fading. 

4. Apparatus for communication with a mobile re- 
ceiver over a fading mobile-satellite channel for digital 
coded voice and data comprising means for trellis cod- 
ing said data, means for block interleaving said trellis 
coded data, means for multiple phase-shift-keyed signal- 
ing said trellis coded and block interleaved data with 
constant amplitude and distinct phase spacing in a signal 
set of phases greater than four, means for quadrature 
modulating the multiple phase-shift-keyed trellis en- 
coded data on an IF carrier, means for inserting two 
pilot tones in the modulated IF carrier, and means for 
modulating the modulated IF carrier on an RF carrier 
for transmission through said mobile-satellite channel, 
and including means at said mobile receiver for syn- 
chronous demodulating of the IF carrier comprising 
means for extracting said two pilot tones, means for 
demodulating said IF carrier using said pilot tones to 
coherently detect said trellis coded data, means for 
block deinterleaving said trellis coded data, and means 
for decoding said deinterleaved trellis coded data, said 
demodulating means including means for recovering 
information for extracted pilot tones for use as side 
information in said deinterleaved trellis decoding 
means. 

5. Apparatus as defined in claim 4 wherein said signal 
set is comprised of eight distinct phases. 
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6. Apparatus as defined in claim 4 including means for 
deriving channel state information from said pilot tones 
at said receiver in the form of a fading channel parame- 
ter, p n , for use as side information, z„, in said means for 
trellis decoding, and said decoding means uses a metric 5 
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of the form m(y„, x„; z„)= - | y„-p„x n | 2 for improving 
the error correction performance of said trellis decod- 
ing means in the presence of channel fading. 

* * * * * 
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